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Abstract 

We compute the next-to-leading order (NLO) perturbative QCD corrections 
to the correlators of nucleon interpolating currents in relativistic nuclear matter. 
The main new result is the calculation of the 0{as) perturbative corrections 
to the coefficient functions of the vector quark condensate in matter. This 
condensate appears in matter due to the violation of Lorentz invariance. The 
NLO perturbative QCD corrections turn out to be large which implies that the 
NLO corrections must be included in a sum rule analysis of the properties of 
both bound nucleons and relativistic nuclear matter. 
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1 Introduction 



The study of bound states in QCD is a difficult problem. With more than 30 years of 
research it is clear that the most promising approach to obtain quantitative informa- 
tion on the properties of hadrons are very likely lattice techniques in particular since 
both computer power and computational methods advanced dramatically since their 
first introduction in the early seventies. Lattice results are now available in many 
hadronic channels and further research is being actively pursued [I]. Nevertheless 
analytical nonlattice techniques can be used to verify at least the consistency of some 
models for hadron description and their predictions. The QCD sum rule analysis is 
based on the operator product expansion (OPE) and serves as a rigorous framework 
for many calculations in the theory of hadrons [21 El HI [5] . QCD sum rules are also 
useful for testing some model dependent approaches P, [3 E] such as the MIT quark- 
bag model [9]. A further important problem is the quantitative description of the 
properties of bound nucleons and relativistic nuclear matter within QCD. 

The proton is the most abundant strongly interacting particle on Earth. It has 
played an important part in particle physics since long ago. In the theory of strong in- 
teractions the proton is a bound state of quarks and gluons with three valence quarks 
fixing its discrete quantum numbers. It has been intensively studied during the last 
fifty years. At present one cannot directly compute the properties of the proton ana- 
lytically from QCD even for an isolated proton. The techniques of the QCD sum rules 
provides a powerful tool for the phenomenological analysis of the proton. The prop- 
erties of baryons have been successfully described within this approach when account 
is taken of the leading vacuum condensates [TOj [TTl [12]. The ever improving accu- 
racy of experimental data requires improvements also in the theoretical description. 
Primarily this means that one has to account for the perturbative QCD corrections 
to the coefficient functions of the OPE for baryonic correlators. The next-to leading 
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order (NLO) perturbative QCD corrections have been calculated for the coefficient 
functions of the OPE for the unit operator and the scalar quark condensate in the 
massless quark limit in fT3]. The perturbative corrections have been found to be 
large. The results for the unit operator were generalized to the massive quark case 
in [H] where again large perturbative corrections were found. The properties of the 
proton in vacuum are well studied in QCD although there is room for improvements 
in the numerical accuracy. 

Protons are traditionally used as targets in accelerator experiments as e.g. in the 
electron scattering on iron at DESY. To analyze the data obtained in these exper- 
iments one needs to know the properties of the protons bound in nuclei, or more 
generally of the nuclear medium. Thus there is a considerable interest in computing 
the parameters of the proton medium. The most obvious reason is that protons are 
part of the nuclei which serve as targets in accelerator experiments. The scatter- 
ing on nuclei is different from the scattering on the proton, and this is important 
for the interpretation of the data. One of the best known examples is the EMC ef- 
fect (e.g. [15]). One of the possible theoretical approaches is to use an effective theory 
where the proton in the medium is treated as an effective particle p!6| [T71 [18] . An- 
other approach is to analyze the properties of nucleon medium within the QCD sum 
rule approach in [19]. In this approach the problem of an accurate determination of 
the parton distributions in nuclei ultimately requires the calculation of perturbative 
corrections to the OPE in the framework of QCD sum rules. 

In the present paper we compute NLO perturbative QCD corrections to the corre- 
lators of baryon interpolating currents in matter. We focus on the calculation of the 
0{as) perturbative corrections to the coefficient functions of the bilinear quark op- 
erators that may lead to the emergence of nonvanishing condensates upon averaging 
over the appropriate physical states. We present new results for the coefficient func- 
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tions of the quark operators in the vector representation (1/2, 1/2) of the Poincare 
group: NLO accuracy is achieved in the expansion in the couphng constant of QCD. 
When averaged over the ground state of matter one obtains nonvanishing values of 
the quark operators in the vector representation (1/2, 1/2) when studying the prop- 
erties of the nucleon medium within the QCD sum rules approach. The presence of 
matter violates Lorentz invariance and thus allows for the appearance of a vector con- 
densate averaged over the matter states. The NLO perturbative QCD corrections to 
the coefficient functions turn out to be large in the MS scheme. This means that one 
must account for the perturbative corrections in applications of sum rules analysis' 
of baryon properties in matter. 

2 Basic expressions for the analysis 

The formulation of the OPE analysis is standard by now. In accordance with the 
QCD sum rule approach, we shall calculate the operator product expansion of two in- 
terpolating currents J{x) which have a nonvanishing overlap with the state of interest. 
The OPE for the quantity 



is performed by means of Wilson's operator product expansion. For the analysis 
of the properties of isolated hadrons in the vacuum one then averages the operator 
product over the ground state of QCD or the physical vacuum to obtain the correlation 
function of two interpolating currents J{x) 



in vacuum. The assumption of the QCD sum approach is that the vacuum expecta- 
tion values of the local operators that appear in the OPE (the so-called condensates) 
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(1) 




(2) 



are nonzero. The calculations are done in perturbative QCD and therefore must be 
performed in the region —q^ > 1 GeV^ where perturbative QCD is valid. However, 
in this region, the effective strong interaction constant as is not very small numeri- 
cally [20] . This forces one to calculate the coefficient functions of the operator product 
expansion in perturbation theory at least up to NLO in order to have sensible results. 
One more reason is of course the general property of perturbation theory that only 
at this order of the perturbative expansion can one reliably fix the renormalization 
group scale /i which determines the numerical values of the coupling constant and 
condensates within the OPE. It turns out that the NLO corrections are rather large 
in many hadronic channels. Even in the case of the correlators of the classical quark- 
antiquark currents the perturbative QCD corrections in the standard MS scheme are 
not small [21]. Thus, the calculations of n(q') should be done at least at NLO in 
to have the precision required by modern applications. 
In the original paper of loffe the current 

r] = eabciulC'j^Ubh^l^dc (3) 

was used to analyze the properties of the proton [TO], u and d are light quark fields 
and C is the charge conjugation matrix with the properties C^yJ^C = 7^ and C = 
—C^^ = —C^ = — By using a Fierz transformation the loffe current t]{x) can 
also be rewritten as a linear combination of the two current operators Oi and O2, 

r]{x) = 2(0i - O2) (4) 

where 

Ol = eabc{UaCdb)'y5Uc, O2 = eabc{u^C-f5db)Uc. (5) 

In fact, the operators Oi and O2 form a complete basis for the lowest dimension 
interpolating currents of the proton with no derivative couplings. We generalize 
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(a) (b) (c) 

Figure 1: LO (a) and NLO propagator-type (b) as well as fish-type (c) corrections of 
the baryonic two-point correlator. 

loffe's current by writing a linear combination of the operators Oi, O2 of the form 

J{x) =0i + t02 (6) 

where t is a mixing parameter. 

The topology needed in the calculation of the LO correlator (see Fig. [I](a)) falls 
into the category of the well-known sunset diagrams. These only contain lines that 
connect two vertices [221 123]. Such a topology also appears in the effective gluon 
low energy correlator for heavy quarks below the production threshold that leads 
to the decay of heavy quarkonia into gluon/photons [211 [25]. Sunset diagrams can 
be calculated very efficiently in configuration space. The NLO perturbative QCD 
corrections are of two types. The first correction is a propagator-type correction 
which is not difficult to compute (cf. Fig. [11(b)). The second correction (Fig. [H^c)) 
comes from the diagrams of the fish type and involves the calculation of an irreducible 
two-loop subdiagram. 

For completeness, and for the convenience of the reader, we present general ex- 
pressions for the corrections in configuration space which can be used for a variety of 
interpolating currents. The efficiency of the configuration space approach has already 
been proven in computing NLO corrections to pentaquark correlators [27] . 
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(a) (b) 



Figure 2: propagator (a) and dipropagator (b) correction with open Dirac indices 

First we list the NLO correction to the fermion propagator (see Fig. [2](a)) in 
configuration space. One has 

^(a:)|NLO = ^(x)Ilo {l - Cp^l {f^WY] , (7) 

where Cp = {N^ — l)/2Nc is the Casimir operator of the color group SU{Nc) {Nc = 3 
for QCD) and 5'(a;)|Lo = 7^a;^Fo(a;^) is the LO fermion propagator where -Fo(x^) is 
defined in the Euclidean domain and is given by 

T{z) is Euler's gamma-function. The space-time dimension is parametrized hj D = 
A — 26 throughout. Written in terms of Fq{x'^) one has 

^(a:)|NLO = Fo(x2)7% {l - Cf^^ [f^WY] ■ (9) 

The renormalization scale nx is the appropriate one for calculations in configuration 
space. This choice avoids the appearance of ln(47r) and 7^; (Euler constant) factors 
in configuration space calculations. The relation of fix and the usual renormalization 
scale /i of the MS-scheme is given by fix = fie""^ /2. Note that the NLO fermion 
propagator is gauge dependent even if the complete calculation is gauge invariant. In 
our calculation we have used diagonal or Feynman gauge. 
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The next quantity needed in our calculation is the NLO correction to the propa- 
gator of a pair of fermions. We call this diquark propagator a dipropagator for short 
and denote it by 5*2 (x). The dipropagator is given in terms of a two- loop amplitude 
with open Dirac indices (see Fig. EJ^b)) which requires a genuine two-loop calculation. 
The result for the dipropagator up to NLO reads [27] 

4:71 e 

+ a3rf^®r3«AMX,)} (10) 

where the coefficients ai, bi and 03 are given by 

11 , 1 11 

ai = -l-—e, bi = -l--e, 03 = -- - 

and where 

rr = ^(7^V-7V7'^)- (11) 

We use the standard notation ® for the direct product of two Dirac or color matrices. 
The generators of the color group algebra t"' appearing in the above expression are 
normalised by the condition tr(t"t^) = 6"'^/2. Eqs. ([7]) and ffTOl) allow one to calculate 
the NLO corrections to n-quark(antiquark) current correlators of any composition 
using purely algebraic algorithms without having to compute any integrals. For ex- 
ample, the form flTUl) has been used in [27| to compute the radiative corrections to 
the pentaquark current correlator. 

The above results need to be renormalized. The renormalization can be done in 
configuration space. To renormalize the single propagator one can use multiplicative 
renormalization. The only ingredient needed is the wave function renormalization 
constant of the fermion. The diagrams involving dipropagators result in mixing of the 
operators under renormalization. Mixing is taken into account through a subtraction 
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of the corresponding vertex divergences generated by the operator that can admix to 
the initial current. The general formula reads 

= + (12) 

where i and j are color indices and ip stands for either the up or down quark fields. 
For definiteness we define our a = i/2[7°,7^]. The resuhs are agam given m 
diagonal or Feynman gauge where we emphasize again that the complete result is 
gauge independent. Note that the part proportional to the product of a-matrices is 
gauge independent. 

Before presenting the results of our calculation we want to remark on the renor- 
malization group properties of the operators Oi_2 defined in Eq. ([5]). They represent a 
complete basis of the operators mixing under renormalization and suffice to perform 
the calculation of the baryonic correlators. As mentioned before the operators O12 
form a basis of operators of lowest dimension for the interpolating currents of the 
nucleon. Since their anomalous dimensions are identical at this order they satisfy the 
same renormalization group evolution. One has 

Note that the numerical value of the anomalous dimension is such that the product 
^m(|u)0i.2(/i) with m(/i) is renormalization group invariant at this order of QCD. At 
NLO the operators mix. The two-loop anomalous dimensions have been computed in 
Ref. [26]. These ingredients allow one to compute the necessary correlator functions. 

3 Correlator including the scalar condensate 

In the OPE one computes the contributions of local operators to the correlator func- 
tion. In case of the vacuum correlators only Lorentz scalars contribute which means 
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that the only nonvanishing condensate is of the form (gg) . The standard vacuum con- 
densate contributions have been calculated before in Refs. [13] including the 0{as) 
corrections. They read 

1 1 dxe''^^T{Jix)J{0)}) = ^n,(g2) + UUq') (13) 

with 
and 

Here Q"^ = — g^, (i/jiIj) = {uu) = {dd). 

In matter a new type of condensate (M|g7^g|M) ^ appears where \M) is the 
ground state of the matter. The quantity (M|g7^g|M) violates Lorentz invariance. 
This is expected since the matter itself fixes a special frame. Therefore one needs to 
account for a new operator in the OPE for the baryonic currents. Including the new 
vector operator one now has 

T{J{x)Jm = C,{x') + C,,{x'){qq} + CI^,Jx){q^,q} (15) 

where the coefficient function Cq^^q{x) of the vector operator {q^y^q} is a four- vector. 
We calculate the coefficient function C^^^g{x) at NLO accuracy by using again con- 
figuration space techniques which were developed in a different setting, namely the 
NLO analysis of pentaquark sum rules [27|. 

First we check on known results using configuration space techniques. We split 
the result for the noncondensate contribution into two parts which refiect the two 
ways of how the Dirac indices have been contracted. One has 

W{x) = -4iV,!(Fo(x^))^l + t^){l + ^((/.^x^)^(l + I)-l)}x^,; 

10 



a b c 

a' b' c' 

(a) (b) 

Figure 3: ordinary (a) and crossover part (b) of the Wick contraction for the operators 
Oi^2 in a symbohc representation. The arrows indicate the direction for the quark 
current within the correlator, the dashed hues are inverted by using transposition and 
charge conjugation. 

= -4iVj(Fo(x2))'(l + t^){l + ^g + ln(/i2x2))}x^:^, (16) 

and 

n^(^) = -Arj(Fo(x^))^l+t)2{l + ^((^^xf Q + + 

= -iVj(Fo(x2))'(l + t)^{l + ^(^ + ln(/.^a:^))}a;2:^ (17) 

for the direct and crossover part of the Wick contraction, respectively, as shown 
symbolically in Fig. [21 The singular parts oc 1/e in each of the contributions cancel 
against counter terms in the course of renormalization performed in Eqs. f|T6|) and ffTT]) . 
Note the different dependence on the mixing parameter t in the two parts of n(a;) = 
n°(x)+n^(x). Note also that our techniques allow for the calculation of all condensate 
corrections but require new modules as indicated in Fig. HI These modules are relevant 
for the calculation of the coefficient functions of the scalar operator {ipip} and the 
vector operator {ip'-y^ip}. 

For the scalar condensate we obtain the results 

Wsix) = 4iVj(Fo(a;2))'(l-t)(l + t)|l + ^(/i^x2)^}x2(^V^)/12, 
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Figure 4: new modules for the correlator corrections including condensates: con- 
densate correction (a) and condensate-propagator correction consisting of two dia- 
grams (b) 

nKx) = iVj(Fo(x2))'(l-t)|l + 3t+^(/i^.x2)'(l + 7t)|a;2(^^)/12. 

where we have assumed {uu) = (dd) = (ipip). Adding up the two contributions gives 

Us{x) = m {Fo{x')y (1 - t) js + 7t + ^ {^ilx'Y (3 + 5t)^x'{^PiJ)/12. (18) 

These results are finite and need not be renormalized. After a Fourier transformation 
the results are in agreement with the results for Umiq"^) in Eq. (fT^ obtained by direct 
integration in momentum space. The corresponding spectral density that appears in 
the integrand of the dispersion representation reads 

Ps{s) = ^^(1 - t) {5 + 7t + ^(3 + 5t)} {^^)/12 (19) 

with s = > 0. The result is proportional to (1 — t) and thus vanishes for t = 1. 
The vanishing of the spectral density at t = 1 is a general property of the correlator 
function related to the chiral structure of the current. Indeed, in the massless limit 
where the chiral symmetry is exact the contribution of the scalar quark condensate 
vanishes in all orders of perturbation theory. 
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(a) (b) (c) (d) (e) 

Figure 5: LO contribution (a) and NLO contributions (b-e) to the correlation function 
including a (scalar or vector) condensate 

4 Correlator including the vector condensate 

We now present our new results for the vector condensate. As mentioned before the 
vector condensate violates Lorentz invariance as a manifestation of the presence of 
matter. It can be calculated by using the same set of diagrams as before (see Fig. [S]). 
For the part proportional to the vector quark operator in the OPE we write 

nv'(x) = (^tx^Avix^) + x^^^Bvix^)) (20) 

where Av{x^) = A^{x^) + ^^(a;^), Bv{x'^) = B^{x^) + ^^(x^) and find the renor- 
malized coefficient functions 

A^y{x') = -8m{Fo{x'))\l + t')[l + ^{ln{f,y) + l^], 

B^x') = -4iVj(Fo(x^))'(l + t2){l + ^(ln(/i^.x2) + l)}, 

A^x') = -2m{Foix^))\l + tr[l + ^(\n{^^y) + l)], 

BUx') = -Nj{Fo{x')y{l + tr[l + ^{\n{^^y)-l)]. (21) 

Note that the curly bracket notation {tp'-f'^ip} refers to an operator before averag- 
ing. We emphasize again that the vacuum expectation value of the vector operator 
vanishes, i.e. (ip'-f'^tlj) = 0, while in matter one has (Mltpj^iplM) ^ 0. 
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In momentum space the correlator function is expanded as 

nv(g) = (iq,Ay{q') + q'l.Bviq')) . (22) 
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For the spectral density one obtains PAyi^) = Pa'{,{.s)+ pa^{s) and PByi^) = Pb^{.s 
Pb^{s) where 

8A^c! . 9^ r «s /15 



and £ = ln(/i^/s). As an example we take loffe's current Eqs. ([3]) and (j4]) which is 
obtained from our general current ([6]) by setting t = — 1 and multiplying by a factor 
of 2. Including the noncondensate contribution and the scalar and vector condensate 
contributions the spectral density is now given by 

(24) 



(47r)4'' I TT V12 J) (47r)2^'^'^^ I 27r J 



37r2 

Our results for the vector condensate confirm the LO results given in [281 ISll- The 
NLO corrections to the vector condensate are new. One can see that they are nu- 
merically large in the MS renormalization scheme at the standard value p = ^/s for 
the renormalization scale. The numerical values of the condensates {M\qq\M) and 
(M|g7^g|M) are nonperturbative parameters of QCD that are built into the sum 
rule analysis. Following [281 [29] we take (M|g7^g|M) = u^^pn where is the four- 
velocity of relativistic nuclear matter and Pat is its density. For the contribution of 
the vector condensate we obtain 



Pn 



i{qu) n + ^ f I + + 2s^ (1 + ^ + 



vr V2 /J I vr V 4 
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(25) 



Canonically QCD sum rules are analyzed at a low scale of the order of 1 GeV. The 
running of the coupling as{Mz) = 0.1176 ± 0.002 [30] to this low scale /i = 1 GeV 
results in as(lGeV)/7r = 0.15 ± 0.1. With this value of the coupling constant the 
NLO correction amounts up to 60% of the leading order result. 

The inclusion of terms proportional to the light quark masses rriu^d do not sub- 
stantially change the quantitative results as the masses of light quarks are small |31j . 
Even for exotic strange matter the results are still valid since the s-quark mass 
is still reasonably small. Two recent 0{aj) QCD sum rule determinations give 
m,(2GeV) = 105±6±7MeV [32] andm,(2GeV) = 92±9MeV [33]. These numbers 
are rather close though somewhat smaller than previous results based on r decay 
data [34J. Even if there is not much hope to detect strange matter on earth, strange 
matter can appear as an intermediate state in the high energy collisions of heavy 
ions. In view of such possible applications the inclusion of strange quark mass cor- 
rections is rather topical. We mention that the contribution of four-quark operators 
are also important [35j. Their contribution can be accounted for in the factorization 
approximation. The result reads 

where Lq = In The accuracy of the factorization approximation for four-quark 

operators has been checked in [36] where the configuration space technique was heavily 
used (see also [37t [38]). 

5 Conclusions 

To summarize, we found an important correction to the correlator of baryon currents 
in media which is needed in the analysis of the properties of relativistic nuclear matter 
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and bound nucleons within the QCD sum rule approach. The correction is given by 
the NLO contribution of QCD perturbation theory expansion to the coefficient func- 
tion of the vector condensate in the OPE of the baryon currents and is not smalL It 
amounts to 60% of the leading order term at a low energy scale relevant to the analy- 
sis of nuclear matter and therefore should be taken into account in phenomenological 
applications. 
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